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基于 GRB051117A 的 能 谱 延 迟 及 其 将 发 特性 的 研究 、 


常 学 钊 ， 彭 朝阳 + ， 王道 周 ， 陈 家 明 


( 云南 师范 大 学 物理 与 电子 信息 学 院 ， 云南 昆明 650002) 


摘要 : ”Margutti 等 人 发 现在 8 MISSAE n LEA UE UOS EORR THES ( 峰值 时 间 、 宽 
度 、 上 升 时 间 、 豪 减 时 间 ) 有 较 强 相关 性 ， 随 着 时 间 的 推移 ， 炮 发 变 宽 ， 谱 延迟 变 大 。 我 们 
fe FH B BOR SUB ERE (DCF) 对 一 个 具有 10 MEAN SHE C33 88 ( GRB051117A ) 进 
行 谱 延迟 研究 ， 发 现在 GRB051117A A , MRS BAIN, PERS. LANA, X 
减 时 间 也 存在 较 强 的 相关 性 Bilag x tpeaks Ws trises taecayo TEMBRA Biz A HS BERT 
间 演 化 , 即 随 着 时 间 的 推移 ， 耀 发 的 宽度 越 宽 ， 谱 延迟 越 大 。 本 工作 将 伽 玛 暴 之 间 存 在 的 谱 
延迟 与 汐 发 特性 参数 之 间 的 关系 拓展 到 了 伽 玛 暴 内 部 ,可 以 帮助 我 们 对 多 炮 发 伽 玛 暴 有 更 深 
的 理解 。 并 且 GRB051117A 内 部 的 耀 发 特性 参数 之 间 存 在 的 相关 关系 与 多 个 伽 玛 暴 之 间 基 
至 是 瞬时 辐射 脉冲 高 度 相 似 ， 这 也 为 X 射线 耀 发 与 瞬时 辐射 脉冲 可 能 是 同一 物理 起 源 提供 
支持 。 

关键 词 : 伽 玛 射线 暴 ; X HAIER ; DCF; WER 


(FS NACE CANS SE) 是 近 几 十 年 来 观测 到 的 最 剧烈 的 天 文 现 象 之 一 ， 伽 玛 暴 由 瞬时 加 
射 和 余辉 两 部 分 组 成 。 对 向 玛 暴 的 能 谱 延 迟 的 研究 近 儿 十 年 也 没有 中 断 过 , 伽 玛 暴 的 能 谱 延 
述 指 的 是 伽 玛 暴 高 低能 光子 到 达 时 间 存 在 差异 , 一 般 来 说 高 能 光子 比 低 能 光子 先 到 达 , 这 就 
FEMS IEREXS, 若是 低能 光子 比 高 能 光子 先 到 达 ， 就 形成 了 负 延 迟 。 研 究 伽 玛 暴 的 能 谱 延迟 
对 伽 玛 暴 的 起 源 和 伽 玛 暴 宇宙 学 研究 具有 重大 的 意义 0 。 通 常 计算 伽 玛 暴 能 谱 延 迟 的 方法 
有 拟 合 光 变 曲线 法 和 交叉 相关 函数 法 (CCF). 

CCF 方法 用 来 估计 伽 玛 暴 的 能 谱 延 迟 是 非常 广泛 的 上 9。Noris Ae SIS AE DR T 
延迟 与 伽 玛 暴 的 光度 存在 负 相 关 的 关系 , Ukwatta 等 人 中 通过 扩大 样本 也 得 到 了 相似 的 结论 ， 
通过 确定 了 谱 延 迟 与 光度 的 负 相关 关系 , 可 以 对 伽 玛 暴 的 红 移 进行 估计 , 这 就 使 得 我 们 可 以 
从 伽 玛 暴 谱 延迟 的 角度 来 对 伽 玛 暴 宇 宙 学 进行 探索 ，Hakklia AP BATSE 数据 库 中 的 
2000 多 个 伽 玛 暴 的 谱 延 迟 进行 了 统计 ， 并 给 出 了 伽 玛 暴 谱 延迟 的 正 负 比例 ; Yi 等 人 中 通过 
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对 长 暴 和 短 暴 的 比较 , 发 现 长 暴 中 的 能 谱 延 迟 与 其 持续 时 间 为 正 相 关 关系 , 而 短 暴 中 并 没有 
发 现 这 个 关系 ， 并 对 比 了 相对 谱 延 迟 R = Tt/AT 的 的 分 布 ， 发 现 长 暴 和 短 暴 的 R 的 分 布 形 状 
没有 差异 但 是 R 的 大 小 存在 差异 ，Roychoudhury 等 人 名 分 析 了 一 个 多 脉冲 的 伽 玛 暴 
CGRB060814) ,发现 这 个 暴 的 第 1、2、4 脉冲 是 正 延迟 而 第 3 个 脉冲 是 负 延 迟 ， 并 且 各 脉 
冲 的 光谱 演化 趋势 相同 。 昌 然 造 成 伽 玛 暴 能 谱 延迟 的 原因 目前 还 不 清楚 , 但 是 伽 玛 暴 谱 的 演 
化 中 和 曲率 效应 外 是 被 大 家 普 裔 认可 的 原因 ， 另 外 也 有 人 认为 伽 玛 暴 谱 延迟 是 谱 演 化 和 曲率 
效应 共同 作用 的 结果 站。 
目前 ， 对 于 伽 玛 暴 谱 延迟 的 研究 ， 主 要 是 对 CGRO/BATSE, SWIFT BAT 和 Fermi 的 数 
据 的 研究 ， 也 就 是 说 主要 集中 在 对 伽 玛 波段 瞬时 辐射 脉冲 的 研究 ， 对 于 X REP TE 
发 的 能 谱 延 迟 研究 是 比较 少 。 而 X 射线 余辉 炮 发 与 瞬时 辐射 脉冲 具有 相同 的 起 源 是 被 广泛 
认可 的 。 Burrowst" 提 出 炮 发 是 在 后 期 由 内 部 激 波 产 生 , 这 需要 GRB 中 央 引 擎 的 重新 激活 ， 
许多 学 者 也 得 到 了 类 似 的 结论 23。Chincarini 等 人 1 利用 SWIFT XRT 数据 库 中 获得 的 113 
DER MATAN X 射线 能 谱 带 能 量 下 耀 发 时 间 特 性 的 演化 。X 射线 炮 发 中 同样 存在 谱 
延迟 的 现象 ，Margutti 等 人 号 对 SWIFT XRT 的 8 个 伽 玛 暴 〈 包 含 9 个 耀 发 ) 进行 了 能 谱 延 
迟 的 研究 ,发 现 耀 发 中 的 谱 延迟 与 溜 发 的 特性 参数 之 间 是 正 相 关 关 系 ,， 并 得 到 了 耀 发 中 谱 延 
迟 与 光度 的 关系 , 这 个 关系 与 瞬时 辐射 脉冲 中 的 关系 相似 , 所 以 得 出 炊 发 与 瞬时 辐射 脉冲 县 
有 同一 起 源 机 制 这 个 结论 ; Sonbas 等 人 59 也 支持 X 射线 焰 发 和 瞬时 辐射 脉冲 具有 共同 起 源 。 
但 是 他 们 并 没有 对 多 兆 发 的 伽 玛 暴 能 谱 延 迟 性 质 做 详细 的 研究 -GRB051117A 是 一 个 具有 复 
杂 光 变 曲 线 结 构 的 多 耀 发 伽 玛 暴 , 在 其 前 2000 秒 的 观测 中 至 少 包 含 10 ASE, Goad 5 AU! 
对 这 个 明亮 的 伽 玛 暴 多 波段 进行 了 研究 , 发 现 其 点 源 强度 和 光谱 的 特性 参数 具有 显著 相关 性 ， 
如 光子 指数 、 硬 度 比 和 拐 折 能 量 ，Chincarini 等 人 中 也 对 这 个 暴 的 炮 发 特性 参数 和 光谱 参数 
之 间 的 关系 进行 了 研究 。 然而 他 们 并 没有 系统 讨论 这 个 暴 的 能 谱 延 迟 性 质 , 所 以 我 们 打算 从 
FETS WEIS ITI FA EWI IE ANF (GRB 051117A )。 

SA M TTT EZ EIB FETE EZ MSE, Ej Marguttii 4A 3x 
中 9 ANE BCP FETE FY ER SAR ES BZ TA AER ARE ER, FF HIE GRB051117A 
内 部 耀 发 的 特性 参数 之 间 的 相关 关系 , 与 瞬时 辐射 脉冲 作 直 观 的 比较 , 使 我 们 可 以 对 这 个 多 
次 发 伽 玛 暴 GRB051117A 有 更 深 的 认识 。 


2 样本 处 理 和 分 析 方 法 


XRT 数据 使 用 HEASOFT 软件 包 v6.26.1 和 相应 的 校准 文件 进行 处 理 。 这 个 暴 的 0.3-10keV 
X 射线 余辉 光 变 曲线 来 自 SwifUXRT Io" 5? https://www.swift.ac.uk/)， 本 网 站 可 以 同时 获 
得 0.3-1.5keV 和 1.5-10keV 能 道 的 数据 ， 当 伽 玛 暴 光源 亮度 大 于 儿 个 count 时 ， 数 据 使 用 
Windowed Timing (WT) 模 式 处 理 ， 当 光源 亮度 降低 时 ， 航 天 器 自动 切换 到 Photon Counting 
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(PC) 模 式 来 观测 伽 玛 暴 ， 为 了 确保 我 们 选择 的 耀 发 是 明亮 的 ， 我 们 选择 GRB 051117A 中 WT 
模式 下 的 光 变 曲线 ，GRB 051117A 光 变 曲线 如 图 1 所 示 。 


Count Rate(0. 3-10keV) 


10 
100 1000 
Time since BAT trigger (s) 


图 1 GRB051117A 光 变 曲线 
Fig.1 The light curve of GRB051117A 

Chincarini “AF Norris BOW AIS TARA 10 HEAL, Norris pe BPO A HE Ae A 
25 HORE ERI ESZO up EFTTA), SEQ TA]. UE {EC ESTA], BREAN OPP PE Norris K 
OPN SU F 
I(y=AAexp | = E fort > t, (1) 

(tt) v 

其 中 4 = exp Qu). w= (1/7), t EXE IE REIR], FUE SCH ELS Norris 这 篇 文 
AM, HERI (HB BON 


tpeak =t; + (T112)? 2: 
耀 发 宽度 为 在 耀 发 开始 和 结束 的 两 个 强度 点 Je 处 的 距离 ， 即 
W = taecay + trise = T2(1 + 4u)" aa 
MEACHAM PR EEA 

_ tdecay 一 trise = —1/2 
k= cm um i e ü an)! í = 


上 升 时 间 tise 和 衰减 时 间 taecay 用 w 和 表示 为 
tdecay,rise x zwü t k) (5) 


SWIFT XRT 的 数据 是 离散 的 ， 所 以 我 们 只 能 使 用 离散 的 交叉 相关 函数 (DCF) 来 估计 
能 谱 延 迟 ，DCEF 与 CCF 具有 相同 的 原理 。 由 于 GRB051117A 光 变 曲线 结构 的 复杂 性 ，DCF 
方法 比拟 合 光 变 曲线 法 更 适用 于 这 个 暴 能 谱 延 迟 的 估计 。 我 们 从 Chincarini 等 人 "文章 中 得 
到 了 GRB051117A 的 耀 发 特性 参数 数据 并 用 DCF 来 估计 每 个 耀 发 的 谱 延 迟 ， 结 果 如 表 1 所 
示 ， DCF xc m ae 


NN d) "T 
i=max(1,1-d) *iVit+d 


其 中 x;，yi 是 两 个 光 变 曲线 第 i 个 时 间 片 的 光子 数 ，N 是 光 变 曲线 时 间 片 的 个 数 ，d 是 y 光 
变 曲 线 相 对 于 x 光 变 曲线 的 偏 移 值 ，DCF 是 关于 d 的 函数 。 我 们 得 到 了 DCF 关于 d 的 变化 
| 线 ， 并 用 高 斯 函数 拟 合 了 这 个 DCF 曲线 。 当 两 个 能 道 的 光 变 曲线 形态 相差 很 大 的 时 候 ， 

我 们 需要 用 更 复杂 的 函数 去 拟 合 DCE 曲线 ， 比 如 高 阶 的 多 项 式 或 者 Norris PAA. ERN 


DCF(d,x,y) = (6) 


T 
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的 样本 中 ， 高 斯 函数 已 经 可 以 较 好 的 拟 合 我 们 的 样本 ， 所 以 我 们 统一 使 用 高 斯 函数 去 拟 合 
DCF 曲线 。 高 斯 拟 合 曲线 的 峰值 是 我 们 所 求 的 两 个 光 变 曲线 相关 性 最 好 的 点 ， 我 们 把 它 设 
为 dd ， 所 以 能 谱 延迟 定义 为 : lag = d' x At, At=T/ANT 是 时 间 序列 的 总 长 度 。 这 种 计算 延 
迟 的 方法 其 实 是 计算 了 整个 耀 发 的 综合 延迟 。 

我 们 用 Monte Carlo 模拟 来 估计 延迟 的 误差 ， 具 体操 作 步 又 如 下 :我 们 假设 光 变 曲线 每 个 
时 间 片 光子 计数 率 的 误差 符合 期 望 为 0 方差 为 1 的 正 态 分 布 , 在 这 个 分 布下 我 们 随机 抽取 每 
个 时 间 片 的 光子 计数 率 ， 产 生 一 组 模拟 光 变 ， 用 DCF 估计 这 组 横 拟 光 变 的 能 谱 延迟 得 到 一 
个 延迟 值 。 重 复 这 个 步骤 1000 次 ， 就 可 以 得 到 1000 个 延迟 值 ， 统 计 这 1000 个 延迟 值 可 以 
得 到 它们 的 标准 偏差 ， 我 们 把 这 个 标准 偏差 当 作 谱 延迟 的 误差 品 。 


3 结果 


3.1 谱 延 迟 与 耀 发 特性 参数 相关 关系 研究 


GRB051117A10 个 耀 发 的 延迟 结果 如 表 1 所 示 , 其 中 , 第 五 个 次 发 在 考虑 误差 的 情况 下 我 
们 认为 延迟 为 0， 其 他 9 个 耀 发 的 谱 延 迟 与 自身 沭 发 特性 参数 的 关系 如 图 2 所 示 。 其 中 ， 红 
色 点 是 GRB051117A 的 9 个 耀 发 ， 蓝 色 三 角 、 黑 色 方块 、 绿 色 星 号 都 是 来 自 Norris 等 人 ”05 
年 的 (以 下 简称 N05) 文章 中 各 能 道 的 脉冲 , NOS 使 用 的 的 是 BATSE/CGRO 的 数据 。 我 们 定义 
lagi = tpeak() 一 tpear()， 其 中 i=2、3、4， 这 样 可 以 计算 出 第 二 〈50-100keV) 、 第 三 

(100-300keV) 、 第 四 〈>300keV) 能 道 与 第 一 〈25-50keV) 能 道 之 间 的 谱 延 迟 ， 因 为 N05 
样本 给 出 的 都 是 长 延迟 暴 (lagu 普遍 >1-2s) ， 所 以 我 们 计算 的 谱 延 迟 主 要 集中 在 10'-10' 
这 个 量 级 之 间 。 之 后 探究 其 与 第 一 能 道 的 脉冲 特性 参数 之 间 的 相关 关系 。 然后 与 漆 发 中 存在 
的 延迟 与 参数 之 间 的 关系 作 比 较 ， 可 以 为 粹 发 与 脉冲 可 能 是 同一 起 源 提供 证 据 。 

图 2 的 左上 部 分 显示 了 GRB051117A 中 谱 延 迟 与 漆 发 峰值 时 间 之 间 的 关系 ， 最 佳 拟 合 关 
系 为 tyeax = 10C64+018)1ag(033+035)， 斯 皮尔 曼 相 关系 数 为 0.45 (p=2.2x10")， 这 说 明 随 着 
时 间 的 推移 ， 炮 发 的 谱 延 迟 越 来 越 大 。 

图 2 的 右上 部 分 显示 了 GRB051117A 中 谱 延 迟 与 峰值 宽度 之 间 的 关系 ， 最 佳 拟 合 关系 为 
w = 104.85+0D1ag(052+02)， 斯 皮尔 曼 相 关系 数 为 0.57 (p=1.2x10")。 在 多 炮 发 伽 玛 暴 中 ， 
泡 发 的 谱 延 迟 的 大 小 与 自身 宽度 成 正比 ， 光 发 越 宽 ， 谱 延迟 越 大 。 

我 们 发 现 ， 这 9 个 耀 发 的 谱 延 迟 中 值 为 2.01s， 平 均值 为 3.55s， 宽 度 的 平均 值 为 130s， 我 
们 之 前 统计 了 48 SARIS CL 137 SAO)» RI 123 个 正 延迟 耀 发 的 谱 延 迟 中 
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s。 粮 发 中 延迟 与 持续 时 间 / 宽 度 之 间 是 


正 相 关 关 系 。GRB051117A 中 的 谱 延 迟 比 较 小 ， 是 因为 它 的 持续 时 间 / 宽 度 比 较 小 。 


图 2 的 左下 部 分 和 右 下 部 分 分 别 显 示 了 炮 发 的 上 升 时 间 和 衰减 时 间 与 谱 延 迟 的 关系 , 可 以 


看 出 , 姿 发 的 上 升 时 间 和 衰减 时 间 都 对 谱 延 迟 的 大 小 产生 影响 。 耀 发 的 上 升 时 间 与 谱 延 迟 的 
相关 关系 为 bise = 1004009) qq 0494018), IE IK SAK ABN 0.8 (p=1.3x107); MERK 
衰减 时 间 与 谱 延 迟 的 相关 关系 为 taecay = 100644013) [qq 0554024) , 斯 皮尔 曼 相 关系 数 为 0.6 
(p=9.68x10”)。 在 GRB051117A 中 炊 发 的 上 升 时 间 与 衰减 时 间 越 大 ， 浪 发 的 谱 延 迟 越 大 ， 


这 与 多 个 伽 玛 暴 之 间 存 在 的 关系 是 一 致 的 "9 。 


从 图 2 可 以 看 出 ， 在 瞬时 辆 射 脉冲 中 ， 也 存在 与 GRB051117A 耀 发 中 类 似 的 延迟 与 参数 之 
间 的 相关 关系 ， 即 在 瞬时 辆 射 脉冲 中 ， 也 普 志 存在 lag X peak W> rises 、taecay 这 样 的 相关 
关系 ， 这 也 为 X 射 线 焰 发 与 瞬时 辐射 脉冲 可 能 是 同一 物理 起 源 提供 支持 。 


3.2 耀 发 特性 参数 之 间 的 相关 关系 研究 


通过 对 GRB051117A 内 10 个 耀 发 的 谱 延 迟 与 特性 参数 之 间 的 关系 的 研究 , 我 们 认为 伽 玛 


其 之 间 存 在 的 延 述 与 特性 参数 之 间 的 关系 可 以 拓展 到 


| 伽 玛 暴 内 部 , 本 节 我 们 将 讨论 多 耀 发 人 


玛 暴 内 部 的 小 发 特性 参数 之 间 的 相关 关系 , 与 多 个 伽 玛 暴 之 间 、 瞬 时 辐射 脉冲 中 存在 的 特性 


ZB BHO REE, MERIT POSER 


上 玛 暴 认识 更 全 面 。 


我 们 从 Margutti 等 人 0512010 年 文章 (以 下 简称 M10) 中 获得 了 9 个 炊 发 的 特性 参数 数据 ， 
从 NOS 中 获得 了 115 个 脉冲 的 特性 参数 数据 ，N05 中 的 脉冲 分 为 4 个 能 道 , 我 们 通过 公式 5 


计算 了 这 些 脉 冲 的 上 升 时 间 tise 与 衰减 时 间 taccay 
具体 定义 请 参考 Norris 等 人 M2005 年 的 工作 。 


图 3 的 左上 部 分 显示 了 粮 发 /脉冲 上 升 时 间 与 衰减 时 间 的 相关 关系 ，GRB051117A 中 次 发 


衰减 时 间 与 上 升 时 间 最 佳 拟 合 关系 为 taecay = 10(036+04Dt@33+026)， 斯 皮尔 曼 相 关系 数 为 


0.75 (p=1.84x10”)， 可 以 看 出 ， 不 管 是 GRB051117A 


、M10 还 是 NOS， 上 升 时 间 都 与 衰减 时 


HEERA, 而 且 所 有 点 都 遵循 tdecay> bise» 在 GRB051117A "m, taecay~ 1.9trise; 前 人 
发 现在 耀 发 中 taccay 大 约 是 2 倍 的 cu, 3503， 这 个 比例 与 本 文 给 出 的 GRB051117A FETE HOHE 


发 衰减 时 间 与 上 升 时 间 的 比例 是 相似 的 。 


图 3 的 右上 部 分 显示 了 峰值 时 间 与 宽度 的 相关 关系 , GRB051117A 中 峰值 时 间 与 宽度 最 佳 


拟 合 关系 为 tyeax = 10(058+05w(0.53+0.19)， 斯 皮尔 曼 相关 系数 为 0.73 (p=2.12x102)， 可 以 


看 出 ，GRB051117A 中 兆 发 峰值 时 间 与 炊 发 宽度 成 


FE 相关 关系 ， 即 随 着 GRBOSIII7A 余辉 


时 间 的 推移 ， 新 产生 的 耀 发 宽度 越 来 越 大 ， 这 可 能 也 是 上 文 提 出 随 着 时 间 的 推移 ， 浪 发 谱 延 


IRIS KIRA. £a 2 以 及 图 3 分析， 在 GRBSII 


17A 中 ， 潍 发 是 随 着 时 间 演 化 的 ， 随 着 


时 间 的 推移 , 炊 发 的 宽度 越 大 , 炊 发 的 延迟 越 大 ,这 一 点 也 支持 了 Margutti 等 人 中 的 结论 。 
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UAH, EP tite ESS Ae Az UR EY Te] Sy SE REEF CRA, ANG GRB051117A 
与 MIO PS ME ACA EK, icr p SCRE Ac WE (LI TR] Sy OE BEA EH ORO AR LCBO FL, TME 
发 中 却 比 较 收敛 ， 这 种 现象 产生 的 原因 是 值得 我 们 以 后 研究 的 问题 。 


图 3 的 左下 和 右 下 部 分 分 别 显示 了 峰值 时 间 与 上 升 时 间 与 衰减 时 间 的 关系 , GRB051117A 


中 峰值 时 间 与 上 升 时 间 最 佳 拟 合 关系 为 tyise = 100155039 OED, 斯 皮尔 曼 相 关系 数 为 
0.78 (p=1.17x102)， 峰 值 时 间 与 衰减 时 间 最 佳 拟 合 关 系 taecoy = 10003550600 72029, gr 
皮尔 曼 相关 系数 为 0.7 (p=3.11x102)。MI10 与 N05 样本 也 存在 与 GRB051117A 中 相似 的 关 


系 o ie ACH EE 时 间 trise 和 衰减 时 间 tdecay 都 与 峰值 时 间 tpeak 有 比较 强 的 相关 性 , 它 们 同 影 


啊 了 峰值 时 间 ， 经 过 计算 ， 


在 GRB051117A 中 ，trise~0.06tpoa tdecay~0.12tpeaxo Chincarini"! 


的 图 5 显示 ， 在 113 ARRP, tue O.05tpeaks tdecay© 0.14tpcar，GRB051117A 也 具有 相似 的 


比例 。 
图 4 显示 了 宽度 与 上 升 、 
最 佳 拟 合 关系 为 w = 100 


173031), 0155105). 斯 皮尔 曼 相 关系 数 为 0.79 p=9.8x103) , HER 


衰减 时 间 之 间 的 相关 关系 , GRB051117A HHA HEE EFT [RI 


宽度 与 泡 发 衰减 时 间 的 最 佳 拟 合 关系 :w = 10 (7039201 107) ,斯 皮尔 曼 相关 系数 为 0.98 
Cp ~ 0)。 宽 度 与 上 升 时 间 有 非常 强 的 相关 性 ， 类 似 的 关系 在 瞬时 辐射 脉冲 中 也 存在 ， 比 如 
Kocevski 等 人 P312003 年 文章 的 图 10 便 给 出 了 脉冲 的 上 升 时 间 与 脉冲 半 宽 具有 非常 强 的 正 相 


关 关 系 。 而 我 们 发 现 不 管 是 伽 玛 暴 之 间 的 耀 发 ， 伽 玛 暴 内 部 的 粮 发 还 是 瞬时 辐射 脉冲 之 间 ， 
都 有 相似 的 宽度 与 上 升 、 衰 减 时 间 之 间 的 相关 关系 。 


图 5 显示 出 的 特性 参数 


(teats Ws Lise laesa) 与 其 不 对 称 性 的 散 点 图 ， k 越 大 的 方向 即 


为 不 对 称 性 越 强 的 方向 。 无论 是 瞬时 辐 册 脉冲 ， 伽 玛 暴 之 间 的 漆 发 和 伽 玛 暴 之 内 的 次 发 我 
们 并 没有 找到 任何 的 相关 性 ， 不 对 称 性 似乎 是 与 其 他 特性 参数 无 关 的 。 


4 总 结 与 讨论 


我 们 使 用 DCF 方法 对 一 


个 具有 10 ATE AZ IT] ESSE CGRBOSIII7AO. 的 谱 延 迟 性 质 


进行 分 析 ， 我 们 发 现在 GRB051117A 中 ， 谱 延迟 的 大 小 与 多发 的 峰值 时 间 呈 正 相 关 关 系 ， 


即 谱 延迟 大 小 随 着 时 间 的 


E 移 呈现 出 整体 上 升 的 趋势 。 有 影响 GRB051117A PRAI EIR 


大 小 的 主要 因素 是 耀 发 本 身 的 宽度 ， 即 耀 发 越 宽 ， 谱 延迟 越 大 。 这 支持 了 Margutti ANI 


WA, BIE AC AE Be a EY 


昌 演 化 的 ， 随 着 时 间 的 推移 ， 焰 发 的 宽度 越 大 ， 谱 延 人 运 也 越 大 。 


Norris 等 人 2005 年 文章 的 


KRR, 宽度 与 脉冲 峰值 时 间 也 呈现 正 相关 关系 , 即 在 瞬时 辐射 脉冲 中 也 存在 随 着 时 间 的 推 


4 指出 号 ， 在 瞬时 辐射 脉冲 谱 延 迟 与 脉冲 宽度 也 呈现 一 个 正 相 
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移 ， 脉 冲 宽度 越 大 ， 谱 延迟 越 大 的 现象 ， 并且 我 们 的 图 2 也 显示 出 脉冲 与 炊 发 具有 相似 的 延 
述 与 特性 参数 之 间 的 相关 关系 ， 这 些 也 为 瞬时 辐射 脉冲 与 X 射线 余辉 耀 发 可 能 具有 相似 产 
生机 制 提供 支持 。 

然而 ， 我 们 发 现 ， 虽 然 X 射线 炮 发 的 宽度 、 上 升 时 间 、 豪 减 时 间 都 大 于 来 自 N05 的 脉冲 ， 
但 是 它们 的 谱 延 迟 却 相差 不 大 , 可 能 造成 这 种 情况 的 原因 有 以 下 两 个 : (1) 样 本 的 选择 效应 ， 
N05 文中 使 用 的 样本 是 大 延迟 脉冲 ， 平 均 延 迟 大 于 1-2s， 而 我 们 使 用 的 GRB051117A WEA 
是 属于 比较 窗 的 泡 发 ， 谱 延迟 比较 小 ; (2) 能 道 选择 不 同 ， 计 算 X 射 线 疙 发 谱 延 迟 的 能 道 是 
相 邻 能 道 ， 而 瞬时 辐射 脉冲 能 道 选择 并 不 是 相 邻 能 道 。 通 常情 况 下 ， 相 邻 能 道 计算 的 谱 延 迟 
结果 比较 小 。 但 根据 图 2 中 显示 的 趋势 , 即使 包括 了 谱 延 迟 较 小 的 暴 ， 瞬 时 辐射 脉冲 也 不 是 
X 射线 耀 发 的 延伸 ， 这 表明 GRB051117A 的 X 射 线 耀 发 与 大 延迟 脉冲 可 能 产生 于 不 同 的 物理 
条 件 ， 不 过 这 也 可 能 只 是 受 样本 选择 的 影响 ， 还 需要 以 后 扩大 样本 来 进行 更 深入 的 研究 。 

GRB051117A 中 净 发 的 谱 延 迟 与 炮 发 上 升 时 间 、 衰 减 时 间 都 呈正 相关 关系 ， 即 耀 发 的 上 
升 时 间 与 衰减 时 间 共 同 影 响 了 谱 延 迟 大 小 。 这 些 结果 与 Margutti 等 人 1 的 图 7 相似 。 
GRB051117A 中 多发 上 升 时 间 与 衰减 时 间 是 一 个 正 相 关 关 系 ，taocay~1.9tise， 这 与 Chincarini 
等 人 MJ 统计 的 113 个 泡 发 之 后 显示 的 结果 相似 。 通 过 研究 GRB051117A 谱 延迟 与 耀 发 特性 
参数 的 关系 ,所 以 我 们 认为 多 个 伽 玛 暴 之 间 存 在 的 谱 延 迟 与 煽 发 特性 参数 之 间 的 关系 也 存在 
FB he SMFS AE Ao 

MER A BADE AC FETE IE ACRE ES ZT) PO AR HE EAL Margutti 等 人 
25 E EAM ELE EEE SZ IRI A DES AR PES ES 58] cp at HE— 
致 ， 这 为 瞬时 辐射 脉冲 与 X 81 2 AS TAS H Bee I6] — YE e EH o 

我 们 估计 的 GRB051117A 中 的 谱 延 迟 相 比较 Margutti 等 人 中 给 出 的 8 个 伽 玛 暴 的 9 RE 
发 的 谱 延 迟 偏 小 ,所 以 导致 了 谱 延 迟 与 宽度 、 峰 值 时 间 、 上 升 时 间 、 衰 减 时 间 相 关 性 虽然 一 
致 ， 但 是 我 们 给 出 的 斜率 却 与 Margutti 等 人 中 图 7 给 出 的 的 斜率 相差 比较 大 ， 造 成 这 种 现 
象 的 原因 可 能 是 : (1) 样本 处 理 的 信 噪 比 与 能 道 选 择 不 同 。 我 们 的 样本 的 信 噪 比 为 2.4， 而 
Margutti 等 人 5 选择 的 信 噪 比 为 4， 这 势必 会 造成 谱 延 迟 结 果 有 偏差 。 我 们 两 个 能 道 为 
0.3-1.5keV 和 1.5$-10keV， 这 是 两 个 相 邻 能 道 ， 由 这 两 个 能 道 计 算 的 谱 延 迟 会 比 Margutti 等 
人 中 文中 0.3-1keV 与 2-3keV 两 个 能 道 计算 的 谱 延 迟 要 小 。(2) 计算 谱 延 迟 的 方法 不 同 , 我 
们 使 用 DCF 方法 估计 谱 延迟 ， 而 Margutti 等 人 号 使 用 了 拟 合 光 变 曲线 比较 te 差 值 的 方法 ， 
这 两 种 方法 计算 出 的 延迟 大 小 应 该 相差 不 大 ， 然 而 ， 对 于 GRB051117A 这 种 具有 极其 复杂 
结构 的 多 炮 发 伽 玛 暴 来 说 ， 在 使 用 Noris 模型 拟 合 炮 发 的 过 程 中 ， 拟 合 优 度 并 不 能 保证 ， 
这 就 造成 了 对 于 tea 的 判定 并 不 准确 ， 导 致 了 谱 延 迟 计算 的 不 准确 ， 这 可 能 也 是 
GRB051117A 某 些 次 发 特性 参数 误差 比较 大 的 原因 。 使 用 DCF 计算 整个 耀 发 的 综合 延迟 会 
避免 这 种 情况 发 生 ， 但 是 由 于 SWIFT XRT 数据 是 离散 的 ，DCF HH ARETE TRA AER, dh 
计 的 准确 性 与 数据 本 身 有 关 , 我 们 使 用 高 斯 模型 拟 合 DCF 曲线 可 以 去 除 噪声 的 干扰 ,然而 ， 
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FE 也 起 着 决定 性 的 因素 ， 比 如 ， 兆 


发 的 上 升 时 间 时 间 间 隅 稳定 在 2s 左右 ， 下 降 时 间 的 时 间 间 隔 突然 变 大 ， 达 到 10s 左右 ， 这 
就 势必 导致 DCF 估计 延迟 不 会 很 准确 。 根 据 不 同 的 样本 ， 选 择 不 同 的 方法 计算 谱 延 迟 ， 这 


样 才 能 使 谱 延 迟 的 研究 更 为 严谨 和 六 


EAA o 


10* 10° 
<7 Lag, of pulses 7 Lag, of pulses 
O Lag,, of pulses D] Lag,, of pulses 7 a * 
3 Lag, of pulses 网 3  Lag,, of pulses a 
10°F] 。 Flares of GRBO50117A * Flares of GRB050117A 
——— Best fit ine g ——— Best fit ine 
S 
C 5 
pus E 
i = Ww z 
10' 
10! 
10° 10° 
10? 107 10? 10! 10? 10? 107 10? 10! 10? 
Lag (s) Lag (s) 
10? 10° 
V Lag,, of pulses * V Lag,, of pulses 
D] Lag, of pulses O Lag, of pulses 局 回 * 
3 Lag, of pulses 米 Lag,, of pulses 
* Flares of GRB050117A * Flares of GRB050117A 
Best fit line 402 Best fit line 
i a 8 
=. 40! $ 
8 
E] EB 
J 3 
10! dien. P 
Sa BH 
eee 
E—SgRCH 
10° 10° 
10? 107 10? 10! 10? 10? 107 10? 10! 10? 
Lag (s) Lag (s) 
BE nt 间 | 特 性 参 尖 263 = BA ft, 3 T. a 品 || EL pie y eu — Au> P i uM. 
图 2 谱 延 迟 与 时 间 特 性 参数 的 相关 关系 。 蓝 色 三 角 、 黑 色 方 块 、 绿 色 星 号 分 别 是 脉冲 第 二 能 道 、 第 三 能 道 、 第 四 能 道 相 比 第 


能 道 的 谱 延 迟 。 红 色 点 是 GRB051117A 的 9 ASHER, 
是 最 佳 拟 合 线 ; tuu, = 10064018) q g(033+035), 


红色 实 线 是 最 佳 


的 相关 关系 ， 纤 
发 衰减 时 间 的 


H 
红色 实 线 是 最 佳 拟 合 线 。 左 上 : 延迟 与 次 发 峰值 时 间 的 相关 关系 ， 


红色 实 线 


斯 皮尔 曼 相关 系数 为 0. 45 (p=2. 2x10")。 右 上 : 延迟 与 六 发 宽度 的 相关 关系 ， 
WER: width = 10(+35+0D1qg(052+02)， 斯 皮尔 曼 相 关系 数 为 0.57 (p=1. 2x10 0. Æ P: HEIR SHER EFTTA 
- 色 实 线 是 最 佳 拟 合 线 : tuu, = 1004+009)1ag049+018)， 斯 皮尔 曼 相 关系 数 为 0.8 (p=1. 3x10 )。 右 下 : BR E 
HRA, SOA ERA ABE: taecay = 10064412) [ag O554024) , WE IK BH RBA 0.6 (p=9. 68X10"). 


Fig.2 The correlation between the spectral lag and the time characteristic parameters. The blue triangle, black square, and green asterisk 


are the spectral lags of the second, third, and fourth channels compared to the first channel of the pulse. The red dots are the ten flares of 


GRB051117A, the solid red line is the best fit line. Upper left: Correlated between lag and t5 , the red line indicates the best-fitting 


relation: tpeak = 1064918) [ag 334035), the spearman correlation coefficient is 0.45 (p=2.2x 10"). Upper right: Correlated between 


lag and width, the red line indicates the best-fitting relation: width = 100850.) [qg(52+02), the spearman correlation coefficient is 0.57 


(p=1.2x10"'). Lower left: Correlated between lag and tye the red line indicates the best-fitting relation: tise = 1004+009) lag 0049+0.18)， 


the spearman correlation coefficient is 0.8 (p21.3x10?). Lower right: Correlated between lag and taecay the red line indicates the 


best-fitting relation: taecay = 100-5520:219 900555929). the spearman correlation coefficient is 0.6 (p-9.68x 107). 
表 1 GRBO51117A10 个 次 发 的 Norris 模型 拟 合 参数 以 及 谱 延 迟 
Table. 1 Norris model fitting parameters and spectral lag of 10 flares of GRBOS1117A 


flare tpeax (8 ) w(s) k trise(S) taccay(S) Lag(s) 
1 145+2.5 97.1+5.4 0.3+0.1 33.1+4.6 64+2.9 4.11+0.26 
2 327.5+- 46+8 0.1+0.2 20.7+3.6 25.4+8.6 0.89+0.07 
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3 3747.8 112461.8 0.6x0.2 22.446.5 89.6457.9 1.48x0.1 
4 437.844.4 58.6414.3 0.2+0.2 24.347.4 34.2412 1.63x0.1 
5 499.1x6.6 85.9429.5 0.4x0.2 2646.8 60428.4 -0.01+0.01 
6 619.6+- 108.2416.7 0.1x0.1 47.5x12.5 60.848.3 1.95x0.11 
7 962.1+4.9 224+18.5 0.2+0.1 90.8+4.8 133.2417.2 9.68x0.68 
8 1104.343.8 109+10.9 0.20.1 43.745.9 65.347.8 2.7840.13 
9 1332.9+2.1 2117.9 0.5x0 56.742 154.3x8.2 7.4340.51 
10 1569+7.3 257.2419 0.40.1 77.749.3 179.5x20.6 2.0140.15 
4 4 T 
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Fig .3 Scatter plot of the temporal parameters. The red dots are the ten flares of GRBO51117A, the black dots are the 9 flares of 


M10, the blue triangle, the black square, the green asterisk, and the black circle are the pulses from the first channel, the second channel, 
the third channel, and the fourth channel in NOS. The solid red line: the best fit line of 10 flares of GRB051117A. Black dotted line: 
tuccay=trise. Upper left: The best fitting relation between the rise time and the decay time of flares: taecay = 10(0365041) 0981026) the 
spearman correlation coefficient is 0.75 (p=1.84 107). Upper right: The best fitting relation between the peak time and the width of flares: 


teak = 10(058+059w(053+019), the spearman correlation coefficient is 0.73 (p=2.12x 107). Lower left: The best fitting relation between 


= (0.15+0.38) ,.(0.5240.17) 
= 10 sn 


the peak time and the rise time of flares: trise , the spearman correlation coefficient is 0.78 (pz1.17x10?). 


540.64) (0552023) 
peak 


Lower right: The best fitting relation between the peak time and the decay time of flares: taecay = 1008 , the spearman 
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correlation coefficient is 0.7 (p=3.11x 10°). 
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Fig 4. Correlation between width and rise/decay time. The red dots are the ten flares of GRB051117A, the black dots are the 9 


flares of M10, the blue triangle, the black square, the green asterisk, and the black circle are the pulses from the first channel, the second 
channel, the third channel, and the fourth channel in NOS. The solid red line: the best fit line of 10 flares of GRBO51117A. Left: The best 
fitting relationship between width and rise time of flares: w = 19 0-1740.31) ¢(0834015)_ the spearman correlation coefficient is 0.79 
(p=9.8x 10"). Right: The best fitting relationship between width and decay time of flares: w = 1060394015) 97) the spearman 


correlation coefficient is 0.98 (p ~ 0). 
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圈 分 别 是 来 自 N05 中 的 第 一 能 道 、 第 二 能 道 、 第 三 能 道 、 第 四 能 道 


到 


Fig.5 Correlation between peak time, width, rise time, decay time and asymmetry k. The red dots are the ten flares of GRB051117A, 
the black dots are the 9 flares of M10, the blue triangle, the black square, the green asterisk, and the black circle are the pulses from the 
first channel, the second channel, the third channel, and the fourth channel in NOS. The solid red line: the best fit line of 10 flares of 


GRBO051117A. 
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Study on spectral lag and flare characters based on 
GRBO51117A 
CHANG Xue-zhao, PENG Zhao-yang, WANG Dao-zhou, CHEN Jia-ming 


(College of Physics and Electronic Information, Yunnan Normal University , Kunming 650092, China ) 


Abstract: Margutti et al. found a strong correlation between the spectral lag of flare and its 
characteristic parameters (peak time, width, rise time, and decay time) in 8 GRBs, and as time 


goes on, the flare becomes wider and the spectral lag becomes larger We use the 
discrete-correlation function (DCF) to study the spectral lag of a multi-flare GRB (GRB051117A ) 


with 10 flares, it is found that within GRBO51117A, the spectral lag is also strongly correlated 


with the peak time, width, rise time and decay time of the flares. this is, 


lag X tpeark, Ws trise» Caecay. The flares will also evolve over time within a GRB. This work 


extends the relationship between the spectral lag between GRBs and the characteristic parameters 
of flares to the interior of GRBs, which can help us have a deeper understanding of multi-flare 
GRBs. Moreover, the correlation between the flare characteristic parameters within GRB051117A 
is highly similar to that between multiple GRBs and even prompt emission pulses. This also 
supports that X-ray flares and prompt emission pulses may have the same physical origin. 

Key words: GRBs; X-ray flare; DCF; spectral lag 


